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The detection of neutrinos from massive stellar collapses can teach us a lot not only about source 
objects but also about microphysics working deep inside them. In this study we discuss quantita- 
tively the possibility to extract information on the properties of dense and hot hadronic matter from 
^— ^ ' neutrino signals coming out of black-hole-forming collapses of non-rotational massive stars. Based 

on our detailed numerical simulations we evaluate the event numbers for SuperKamiokande with 
neutrino oscillations being fully taken into account. We demonstrate that the event numbers from a 
Galactic event are large enough not only to detect it but also to distinguish one hadronic equation of 
state from another by our statistical method assuming the same progenitor model and non-rotation. 
This means that the massive stellar collapse can be a unique probe into hadron physics and will be 
' a promising target of the nascent neutrino astronomy. 

PACS numbers: 26.50.+X, 95.85.Ry, 97.60.-s, 21.65.Mn 
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£!(■ I- INTRODUCTION 

© ■ 

One of the important roles of astrophysics is to explore physics under extreme conditions that are difficult to realize 
in terrestrial experiments. In this sense, hadron physics at supra- nuclear densities and non-zero temperatures (e.g., 
hyperon appearance, quark deconfmement and so on) is a natural target of astrophysics and the gravitational collapse 
of massive stars at the end of their lives will set the stage [H, Hj]. In particular, black- hole- forming collapses expected 
for very massive stars with a mass larger than ~30 solar masses (Mq) will be the most promising site. Although such 
an event has not been observed yet, a black hole candidate with an estimated mass of 24-33M Q was discovered [3j and 
this might be a remnant of the collapse of such a massive star. Recently, a regular monitoring of ^10 6 supergiants 
within a distance of 10 Mpc is proposed to detect their silent disappearances [4| and the black-hole-forming collapse, 
which would be invisible optically, might be observed that way. 

In our previous studies [5-7], we showed that the event, which we refer to as the "failed supernova" hereafter, 
is as bright in neutrino emissions as ordinary core-collapse supernovae. We showed also that its time evolutions of 
luminosities and spectra are qualitatively different from those of the ordinary supernova explosion and the ensuing 
proto- neutron star cooling Q, which may lead to the delayed black hole formation for some reasons iEl. Our 
numerical data were adopted as a reliable basis to predict the relic neutrino background from stellar collapses (llj . 
More importantly, however, we also demonstrated by employing different hadronic equations of state (EOS) that the 
duration of neutrino emissions from the failed supernova is sensitive to the stiffness of EOS at supra-nuclear densities 
and, therefore, that the observation of neutrinos from such an event will provide us with valuable information on the 
properties of dense and hot hadronic matter as well as on the maximum mass of proto-neutron stars. 

Although this approach is simple and robust, valid irrespective of neutrino oscillations, it can not distinguish EOS's 
with a similar duration of neutrino emissions: a soft nucleonic EOS and a hyperonic EOS, for example. In this study, 
we attempt to break this degeneracy by analyzing more in detail the time variation of neutrino numbers observed at 
a terrestrial detector, which we refer to as the "light curve" hereafter. While we have studied the detection of failed 
supernova neutrinos taking fully into account the neutrino oscillation and its parameter dependence so far we 
innovate a new method here by employing the Kolmogolov-Smirnov (KS) test, which is free from the ambiguity of 
the distance to the progenitor. We adopt the results of our detailed numerical simulations and evaluate the neutrino 
event number for SuperKamiokande (SK) as a representative of currently operating neutrino detectors. This is the 
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FIG. 1: Luminosities and average energies of z/ e for four different EOS's, LSf80-EOS (dash-dotted), LS220-EOS (dashed), 
Shen-EOS (dotted) and Hyperon-EOS (solid). Time is measured from the bounce. 



first ever serious self-contained attempt to demonstrate that for Galactic events it is indeed possible to break the 
degeneracy for hadronic EOS's by the statistical analysis. 

We arrange this paper as follows. A brief review of the neutrino detection and a description of the newly proposed 
statistical method are given in Sec. [XT] The main results of our study are reported in Sec. IIII1 In Sec. IIV1 we mention 
the possible uncertainties and observational issues. Sec. [V] is devoted to a summary. 



II. METHODS 



A. Setups for failed supernova neutrino detection 



The evaluation of the light curve of neutrinos from the failed supernova can be roughly divided into three parts. 
The first step is a computation of the neutrino luminosity and spectrum at the source. The general relativistic 
^-radiation-hydrodynamics code, which solves the Boltzmann equations for neutrinos together with the Lagrangian 
hydrodynamics under spherical symmetry, is utilized to quantitatively compute the dynamics as well as the neutrino 
luminosities and spectra up to the black hole formation. This code passed a couple of well known standard tests and 
a detailed comparison with Monte Carlo simulations fl3T - [l5j ]. The numerical errors are estimated to be ~10% from 
the computations with lower resolutions in Ref. [Hj]. The progenitor model with 4OM [ijj is adopted as the initial 
condition for the dynamical simulations. 

A hadronic EOS is needed at this stage. It should be emphasized that it is not the intention of the paper to endorse 
a particular EOS but that EOS's which are available for astrophysical numerical simulations, that is, subroutines 
or tables that provide thermodynamic variables in wide ranges of density, temperature and proton fraction are very 
limited at present. For example, an EOS table including hyperons has been provided only by Ishizuka et al. flSl so 
far, based on the relativistic mean field theory. The EOS's by Lattimer & Swesty [Ig| and by Shen et al. 
are the ones from the limited options for the nucleonic EOS. The former is an extension of the compressible liquid 
drop model with three choices of incompressibility (K — 180, 220, 375 MeV) and the latter is based on the same 
frame work as in Ref. [l8[ but without hyperons. In this paper, we employ the EOS's by Ref. [l8| (Hyperon-EOS), 
Ref. [13[2l| (Shen-EOS) and Ref. [H with K = 180 MeV (LS180-EOS) and 220 MeV (LS220-EOS). Whereas the 
results for Hyperon-EOS, Shen-EOS and LS180-EOS have already been given also in Ref. [l|, LS220-EOS model is 
newly computed. 

In Fig.Q] we show the luminosities and average energies of electron-type anti-neutrino for these EOS's. The result 
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TABLE I: Event numbers for the progenitor at the distance of 10 kpc. iV all pc and 7V 5 s pc denote the event numbers until the 
end of simulation and up to 0.5 sec after bounce, respectively. 
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for Hyperon-EOS is almost the same as that for Shen-EOS until hyperons appear around ~0.5 sec, since the two EOS's 
are identical in the low density regime. The Shen-EOS, which is the stillest (K — 281 MeV) among these EOS's, is 
sufficiently distinguishable from the other three just by its longer duration (1.345 sec) of neutrino emissions. This 
will not be the case, on the other hand, for Hyperon-EOS and LS-EOS's with K ~ 200 MeV because the duration for 
Hyperon-EOS (0.682 sec) is not very different from those for LS180-EOS (0.566 sec) and LS220-EOS (0.784 sec). This 
is the degeneracy problem mentioned earlier. It is the main purpose of the study to demonstrate that the relatively 
small difference displayed in Fig. [1] is still sufficient to distinguish one EOS from the others by the statistical analysis 
given below. 

The second step for the evaluation of the light curve is to incorporate the neutrino oscillations, which take place as 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect during the propagations through the stellar envelope and the Earth. 
Since the reactions at the detector depend on the neutrino flavor, a proper treatment of the neutrino oscillations is 
mandatory. In this paper, we ignore the Earth effect, since the event number will not be affected very much [l^. As 
for the undetermined parameters in the neutrino oscillation, namely the mixing angle, #13, and the mass hierarchy, 
we choose two limiting cases: the normal mass hierarchy with sin 2 #13 = 10 -8 and the inverted mass hierarchy with 
sin 2 #13 = 10~ 2 . For the detailed dependences on the choice of these parameters, we refer readers to Ref. [l2l j. 

The last part of the procedure is an evaluation of event numbers at the detector. In this paper, we take SK from 
currently operating neutrino detectors. The neutrino reactions in the detector that we take into account are as follows: 

v e +p — ► e+ + n, (1) 

u + e — ► v + e, (2) 

v e + 16 — > e+ 16 F, (3) 

v e + 16 — > e+ + 16 N. (4) 

The first reaction gives a dominant contribution to the event number and we take its cross section from Ref. 22j. 
The second reaction occurs for all flavors of neutrino but with different cross sections, which are taken from Ref. 23j . 
The cross sections for the others are adopted from Ref. (24|. We assume that the fiducial volume is 22.5 kton and 
the trigger efficiency is 100% at 4.5 MeV and 0% at 2.9 MeV, which are the values at the end of SuperKamiokande 
I [25|. The energy resolution was 14.2% for E e = 10 MeV at that time [25| and roughly proportional to \/W e 26], 
where E e is the kinetic energy of scattered electrons and positrons. We choose the width of energy bin to be 1 MeV 
in this study. The event numbers for the progenitor at the distance of 10 kpc are listed in Table U for eight different 
combinations of EOS and mixing parameters. 



B. Statistical analysis 



In order to see if two different hadronic EOS's can be distinguished by the neutrino observations, we take the 
following strategy. We first generate the "observational data" by Monte Carlo simulations based on the light curve 
obtained above for one EOS. We then take the light curve for another EOS as the "theoretical" model and employ 
the KS test to judge if the difference is statistically significant. In this study, we adopt Hyperon-EOS for the former 
and LS180- and LS220-EOS's for the latter. 

In our analysis, we utilize the normalized time profiles of the cumulative event numbers for the first 0.5 sec of 
detections so that we could ignore the neutrino emissions after the black hole formation and the uncertainties of the 
distance measurement to the progenitor. In this study, we consider two cases, in which the total event numbers up to 
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0.5 sec, -/V0.5 s , are 10000 and 400. The former roughly corresponds to the source located at the Galactic center while 
the latter represents an event in the Large Magellanic Cloud (LMC) or Small Magellanic Cloud (SMC). 

In the KS test, the so-called KS measure, £>ks, is defined as the maximum difference [13] between the two data: 

D KS = max |/theor.(*) ~ /obs.(*)|) ( 5 ) 
i<0.5s 

where /theor.(i) and / bs.(i) are the theoretical and observational time profiles of neutrino events, respectively. As 
mentioned already, they are normalized as /theor.(0.5s) = Jobs. (0.5 s) = 1 in our analysis. If Dks > 0.01622 (0.0811) 
for 7Vo.5s = 10000 (400), the theoretical model is rejected at the confidence level of 99%. 

We also employ time-shifted observational data, since we can not know the onset of neutrino detections precisely. 
In this case, the KS measure becomes a function of the time shift, i s hift, as 

Acs(*shift) = max |/theor.(t) - fobs.(t - Wft)|, (6) 
t<0.5s 

where the time profiles are normalized as /theor.(0.5 s) = .fobs. (0.5 s — £ s hift) = 1- If the minimum value of the 

KS measure with respect to the time shift satisfies the above criteria, i.e., min DKs(ishift) > 0.01622 (0.0811) for 

*ahift 

-Wo. 5 s = 10000 (400), we can conclude at 99% C.L. that the two data are distinct from each other. 

III. RESULTS 
A. Analysis without time-shift 

We first show the results of the analysis without the time-shift. In Fig. [2J the un-normalized time profiles of 
cumulative event numbers are compared between a MC data ("observational data") randomly picked up from 100000 
realizations for Hyperon-EOS and the "theoretical" estimations for LS180-EOS and LS220-EOS. Since the growth 
rates of neutrino luminosity and mean energy are larger for LS-EOS's than for Hyperon-EOS (see Fig.rjJ, the lines for 
the former are lower than for the latter. The vertical bars are corresponding to ±0.01622 (0.0811) in the normalized 
profiles for Nq^ s = 10000 (400) and the theoretical models outside these error bars are rejected by the KS test at 99% 
C.L. It is clear from the figure that for Nq^ s = 10000 or the Galactic source, we can easily distinguish Hyperon-EOS 
from LS-EOS's. In fact, this is true not only for this particular data but also for all 100000 MC realizations in each 
set of the mixing parameters. Thus we can conclude that Hyperon-EOS and LS-EOS's employed in this study are 
distinguishable for Galactic events. 

The claim depends on the mixing parameters for events in LMC or SMC, that is, for iVo.5 S = 400 as can be 
inferred from Fig.[5J In fact, for the normal mass hierarchy with sin 2 #13 = 10~ 8 , Hyperon-EOS is not distinguishable 
from LS220-EOS (LS180-EOS) for more than 65% (40%) of 100000 MC realizations. In the case of the inverted 
mass hierarchy with sin 2 $13 = 10 -2 , however, the distinction fails 5218 (512) times among 100000 MC realizations 
for LS220-EOS (LS180-EOS). The reason why the latter case of the mixing parameters is easier is as follows. The 
EOS dependence is stronger for muon-type and tau-type neutrinos and their anti-neutrinos than for electron-type 
counterparts (see Fig. 3 of Ref. |1J). A portion of and v r is converted to v e , the dominant contributor to the event 
number via reaction ([T]). While 16% of and v T become v e for the normal mass hierarchy with sin 2 #13 = 10 -8 , a 
half of and D T are converted to D e for the inverted mass hierarchy with sin 2 #13 = 10 -2 , thus yielding the observed 
EOS dependence of the event numbers. 

B. Analysis with time-shift 

Next, we demonstrate that the conclusion is unchanged if the time-shift is included in the analysis. In Fig. [3J the 
distributions of the KS measures for 100000 MC realizations are shown as a function of the time-shift, i s hift- LS180- 
EOS (LS220-EOS) is employed to obtain the "theoretical" model in the left (right) panel. For Galactic sources, that 
is, in the case of iVo.5 S = 10000, all of 100000 MC realizations have Dks larger than 0.01622, the lower limit for the 
rejection by the KS test, for any value of the time-shift in the case of the normal mass hierarchy with sin 2 #13 = 10 -8 . 
This is all the more true for the inverted mass hierarchy with sin 2 #13 = 10~ 2 . 

For Nq.5s = 400, however, the above conclusion is somewhat compromised. In the case of the normal mass hierarchy 
with sin 2 #13 = 10 -8 , for example, Hyperon-EOS and LS180-EOS are not distinguishable in more than 90% of 100000 
MC realizations when i s hift = 0.06 sec is assumed. For the inverted mass hierarchy with sin 2 #13 = 10 -2 , on the 
other hand, more than 75% (90%) of MC realizations have a KS measure larger than the critical value for LS220-EOS 
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FIG. 2: Un-normalized time profiles of cumulative event numbers for two mixing parameter sets and two total event numbers. 
The observational data are shown by solid lines whereas the theoretical estimations are displayed by dash-dotted lines for 
LS180-EOS and dashed lines for LS220-EOS. The error bars are corresponding to 99% C.L. in the KS test. 




FIG. 3: The KS measures, -Dks, as a function of the time-shift, tahitt, for LS180-EOS (left panel) and LS220-EOS (right panel). 
50%, 90% and 99% of -Dks's lie above the solid, dashed and dot-dashed lines in 100000 MC realizations, respectively. The 
horizontal dotted lines represent the lower limit to reject the model by KS test. 
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(LS180-EOS). In conclusion, for the event in LMC or SMC, the successful distinction of Hyperon- and LS-EOS's 
depends on the neutrino mixing parameters; the inverted mass hierarchy with sin 2 (9 13 = 10~ 2 is more advantageous. 

IV. DISCUSSION 

Certainly, there are some issues remaining unaddressed. First of all, investigating other EOS's constructed by 
different manners is important future work. In particular, hyperon mixing in neutron star matter has been studied 
by the microscopic approach [28 30] and it is desired to be applied to EOS tables including finite temperature. If 
neutrinos have masses of ~ eV, the arrival time depends on the neutrino energy and the difference becomes ~0.01 sec 
between the neutrinos with an energy of E ~ 10 MeV and those with >30 MeV [31] for Galactic sources. Recently, 
the core collapse simulations for the same progenitor model are done utilizing a spherically symmetric ^-radiation- 
hydrodynamics code by another grope [32[ . While they conclude that their results are in qualitative agreement with 
our study, there are some quantitative differences. For instance, in the case with LS180-EOS, they found that the 
interval time between the bounce and black hole formation is 0.436 sec, which is ~20% shorter than our result. Further 
numerical investigation is important. 

Multi-dimensional hydrodynamical effects such as rotations and magnetic fields, which are not taken into account 
in our model, may also affect the signals. In the context of supernovae explosions, Ref. [33| compared the neutrino 
signal by an one-dimensional simulation with those by two-dimensional simulations (with and without rotation) for 
the collapse of a 15M Q star adopting LS180-EOS (see Figs. 1 and 14 of this paper). According to this paper, when 
the initial rotation of progenitor core is a constant angular frequency of 0.5 rad s _1 , the stellar rotation does not 
affect the luminosity and mean energy of the emitted neutrinos very much. Incidentally, the neutrino emission after 
the collapse to a black hole is estimated to be negligible for low-angular momentum cases [34] ]. 

Uncertainties with progenitor structures are another concern. While a certain class of progenitors provides similarly 
short neutrino burst, it has been shown that the density profile of outer layer may affect the duration Q. Especially, 
in the case with lower matter density outside the core [35j, the interval time between the bounce and black hole 
formation becomes longer as 1.477 sec even for the model with 40M Q and LS180-EOS 32]. However, if a neutrino 
event is actually detected, we will be able to determine the direction of the progenitor to some extent by the neutrino 
detection itself [36[ and the progenitor is highly likely to be identified or at least constrained by earlier records of 
optical observations as in the case of SN1987A. Then we can study the progenitor dependence much more efficiently. 

As discussed in Refs. [H[l2j], the event rate of the black- hole-forming failed supernovae is estimated to be somewhat 
low ^0.008/yr. This problem may be circumvented, however, by deploying a large detector such as Deep-TITAND 
with a fiducial volume of 5 Mton [37| proposed currently. With this large facility, the event rate goes up to ~0.02/yr 
since we should be able to detect ^400 neutrinos for the first 0.5 sec of a failed supernova in galaxies as far away as 
the Andromeda galaxy (M31) at 780 kpc from us. This is large enough to distinguish the EOS's we studied in this 
letter for the inverted mass hierarchy with sin 2 9-\ 3 = 10~ 2 . Note that the cumulative core-collapse supernova rate 
within this range is 10% of that within 10 Mpc [38|, where a regular monitoring of ~10 6 supergiants is proposed in 
Ref. [|. 

V. SUMMARY 

This study is the first ever serious self-contained attempt to assess the neutrino signals from black-hole-forming failed 
supernovae, which would be observed by the currently operating terrestrial detector. Based on our detailed numerical 
simulations, we have evaluated the event number of neutrinos emitted from black-hole-forming failed supernova for 
some EOS's of nuclear matter. Ambiguities on the neutrino mixing parameters and the onset of the neutrino emission 
have also been taken into account for the evaluation. Assuming the same progenitor model and non-rotation, we 
have shown that we will be able to constrain the EOS of nuclear matter not only from the duration time of neutrino 
emission but also from the time variations of neutrino event number for the progenitor in our Galaxy. Moreover, in the 
case of the inverted mass hierarchy with sin 2 #13 = 10~ 2 , the constraint is favorable even for the progenitor in LMC 
or SMC. The positive results presented here clearly indicate the importance of further investigations of the hadronic 
EOS at supra-nuclear densities based on better formulations and encourage in particular those engaged in the study 
of the EOS of hadronic matter to prepare their latest results in a form available for astrophysical simulations. We 
hope this paper will advance such collaborations further. 
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